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ABSTJIACT 


A pilot study to detect relationships between anomalies in the 
stratosphere and later anomalies at the surface at the same station 
has identified tvjo potential relationships. The anomaly of October 
precipitation at mid- and high-latitude stations appears related to 
the date the fall reversal of zonal v;ind occurs; and the magnitude of 
the anomaly of surface temperature five or six months following a key 
month appears related to the stratospheric zonal wind speed during the 
key month. Both relations are statistically significant at the five 
percent level, and both should be analyzed in more detail for possible 
use in long range forecasting. Horizontal fields from maps should be 
used to detcrm.inc the interlcvel relations on a hemispheric or global 
basis to understand why these relationships exist. 
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I. INTRODUCTION 

The purpose of this study is to identify anoiralous parameters and 
events in the stratosphere which are precursors of anotnalous tropospheric 
seasonal climatic regimes. The present effort is a pilot study to delin- 
eate those stratospheric parameters which are coupled with a later 
tropospheric condition, and to formulate plans for future studies which 
will lead to an understanding of how the stratospheric anomalies develop 
and how they can be used as input for long-range forecasting. 

The importance of reliable long-range forecasts is self-evident; 
however, present numerical weather prediction models are unable to fore- 
cast reliably for more than a fraction of a month. Laurman (1975) has 
reviewed the prospects for extending the forecast period and concludes 
that empirical approaches to long-range forecasting should be the Tiiost 
useful at this time. The empirical approach follcjwed licre is based on 
the fact that the layers of the atmosphere are intimately coupled and that 
events in one layer can be associated with events in another layer. 

Simultaneous, compensating temperature changes in the troposphere 
and stratosphere have been recognized for over fifty years, aud 
a simultaneous relation between the decrease of wind at 45-60 kra at White 
Sands and an increase in surface precipitation over New Mexico has been 
stated by Jetton (1968). Events in the stratosphere which precede events 
in the troposphere are generally less v;ell documented or accepted than 
simultaneous events, but some notable cases are available. For example, 
Drogaytsev (1972) has carefully analyzed the relationship between develop- 
ment of the fall reversal of stratospheric wind and the temperature regime 
at the surface the following winter. Labitzke (1965) has shown that 
sudden warmings in the stratosphere are followed by a blocking circulation 
at sea level in about ten days. 

Identification of a relation between a stratospheric event and a 
later tropospheric event is important, but before the relation can be 
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used to nake forecasts with confidence it must be understood. For 
example, the modeling efforts of Trcnberth (1973) and others have led to 
a relatively good understanding of the processes during a sudden x;arming, 
and so the troposi'iuric blocking after a sudden warning is seen to be a 
natural outcome of the chain of events. Trcnberth states that "blocking 
and the warming may be regarded as part of the sam.e process, not that the 
blocking is som.ehow caused by the warming.' Thus, while a description of 
inter-level relations found in this study is the central theme below, a 
discussion of possible future efforts to gain better understanding of 
those relations is also given. 

II. DATA 

I 

The present pilot study is based on single station analysis of 
stratospheric and tropospheric anomalies. The seven rocket (MR’.s) stations 
listed in Table 1 were selected on the basis of geographical coverage and 
quantity of data available. While lUllI data are available since 1960 at 
some stations, the number and quality of observations increased steadily 
through the early 1960*s and has remained fairly constant to the present. 
Thus, the ten year period 1964-1973 was chosen for the present analysis and 
all anomalies (surface and upper-air) are defined with respect to this 
period. Monthly means of surface temperature and precipitation v:cre taken 
from "Monthly Climatic Data fer the World." Daily radiosonde data, obtained 
from the World Data Center, Asheville, were available through 1970; data 
for 1971-1973 were not received in time to be included in this study. 

The MRN data were generally used at 2 km intervals, 20-60 km, although 
when several adjacent layers indicated a similar relationship with a 
surface parameter the analysis was repeated using average values of variable 
height layers. Monthly means and deviations of temperature, zonal and 
meridional wind speed (T,l',V) \.ere used, along with monthly mean values of 
the vertical shear or lapse rate. Also used were llie dates of the autumn 
reversal of the zonal wind direction, ll^ose dates are given on a zonal 
mean basis in Belmont, ct al. (1975), but for the present study dates at 


the individual stations were required, with the finest tiir.e resolution 
possible. The fall reversal is defined here as the last date the zonal 
wind bocame westerly after having been easterly for at least five consecu- 
tive days. Time-hiight sections of individual daily soundings were plotted 
and the reversal progression analyzed. In most cases, the reversal date 
could be determined within + three days, and in some cases within + one 
day. As discussed in Belmont, et al. (1975), the low-latitude reversal 
is weak and difficult to define. Thus, the reversal analysis was limited 
to stations north of 35°N. 

The spring reversal process is much more erratic than that in the 
fall. Because of nid-wintcr stratospheric warmings the zonal flow may 
experience several changes of direction between January and June. As a 
result of this behavior, "objective" definitions can be made such that 
the spring reversal "date" is forced to correlate with nearly any 
parameter. For this reason the effect of the spring reversal date was 
not examined here. 

III. RESULTS 

A . Relation of the Autumn Reversal of Stratospheric Zonal Wind to 

Precipitation Anoi’.alies 

The anomalies of the autumn reversal dates, relative to the mean 
reversal date at each station, were compared with surface parameters for 
the autumn and winter months which follovjcd. The reversal date anomaly at 
30 km was found to correlate well with the precipitation anomaly in 
October at the three high latitude stations where reversals could be 
unambiguously determined for several years. Note in Table 2 that the signs 
of the anomalies are the same in eight out of ten cases at Cliurchill and 
seven of nine cases at Grccly and Wallops. The chances of this happening 
are 5% and 97», respectively. If the three data sets are assumed indepen- 
dent and combined, the chance of getting twenty-two identical signs out 
of twenty-eight cases is m'ar 1/i. 
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The magnitude of the precipitation anomaly is marginally 
correlated with the fall reversal anomaly for those years when the 
algebraic signs are identical. The linear correlation coefficients and 
their levels of statistical significance are: Greely, 0.71 (57,); 

Churchill, 0.53 (~!07,); V.’allops, 0.67 (5-107.); and all stations, 0.46 (57,) 

In an effort to determine if the above relationship could be 
defined with greater time resolution the weekly precipitation values for 
October at Wallops and Greely were extracted from Monthly Weather Reviev; 
when available. Ko pattern was suggested by this exercise. 

A physical reason for this relation is not easily seen. However 
it is interesting to note that the relation of October precipitation with 
reversal date is much stronger than the correspond: ..g relation of precip- 
itation with the 2C0 mb height or its standard deviation (Table 2). WTiile 
analysis of horizontal fields is necessary for confirmation, this suggests 
that the precipitation anomalies may not be due to a particular standing 
wave-pattern or to enhanced traveling \;aves, but probably to a comhiiiation 
of the phase and amplitude of the predominant standing waves and the 
number of traveling waves. 

The autumn reversal is basically due to cooling near the pole 
and the associated changes in thermal and height gradients. However, a 
cursory review of high-altitude synoptic maps shows that circulation 
features influence the reversal date at a given location. Note that 
Drogaytsev (1972) found the strongest relation between stratospheric 
reversals and the following season’s surface temperature when he strati- 
fied the data according to the location of the earliest reversal. Also, 
he reports that the location of the earliest reversal coincides with the 
location of minimum total ozone. This may be pertinent for unders tand ing 
the present results because of the strong links between the tropopausc 
heiglit and character with total ozone and \7ith weather systems in the 
lower troposphere. Perhaps the synoptic situation in the stratosjvhcre 
near reversal time and the morphology of the tropopausc should ce used to 
further stratify the present study to find a stronger rclationsliip and 
to gain understanding of why the above relation exists. 



Maps based on radiosonde and MRN data are available for the 
entire period used here, and could be analyzed to determine synoptic 
conditions near each fall reversal date. Zonal haromonic analysis can be 
used to find the amplitude and phase of the predominant wavenumbers in 
the stratosphere at reversal time and to tiace these v;aves dovmward into 
the troposphere and fon/ard in time to October. To study several latitudes 
without resorting to spherical analysis, it is possible to stratify the 
years according to curvature of the streamlines at each MRN station at 
the time of fall reversal (i.e., determine if the station is south of a 
low, or north of a high). Lastly, the traveling wave energy and the 
ratio of traveling to standing wave energy should be determined. Because 
the fall reversal is' basically a reordering effect due to radiation losses 
near the pole, its haste (delay) is probably influenced by dynamic cooling 
(heating) due to eddy heat flux. In a similar manner, abnormally large 
amounts of transient v;ave energy in the troposphere, coupled with the 
phase of the predominant standing v;ave, will determine \;hich longitudes 
are apt to receive the greatest precipitation amounts. 

B. Relation of Stratospheric Zonal Wind Sneed to Surface Temnerature 

Anomalies 

One method of identifying potential relationships is the super- 
posed epoch analysis. In the present study a physically significant key 
date (or key month) could not be defined, so the follov;ing procedure was 
used; Each of the twelve months, for the ten years, v;as sorted into tv:o 
classes according to the sign of the anomaly of the monthly mean zonal 
wind speed. Thus, for example, of the ten Januaries, those months whose 
mean wind was higher than the ten year mean, went into the positive group 
for January. This makes twenty-four groups of months. For each of twelve 
calendar months, the ten year mean surface temperature was prepared, and 
the difference, or anomaly, was computed for each of the ten corresponding 
months. The mean of the absolute values in eacli group of positive (and of 
negative) anomalies was next computed, and for each month afterward, for 
the same calendar year sequenro, up to eleven months. The difference of 
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the mean absolute values of the anomalies of surface temperature for each 
of the two groups of months was then determined for each month, and for 
each month af terv^’axds . The results were arranged in a 12 x 12 difference 
array where there are twelve columns for the twelve months, and the rows 
represent months after the starting month. Note that the year groupings 
will be different for each starting month. 

This procedure was repeated for each height interval and at 
each station. The largest value was selected from each difference array, 
and its starting month and number of months after the starting month were 
noted (i.e,, its coordinates in the array). It was found that for the 
height interval 35-45 km the largest value always occurred five or six 
months after the starting month and that the starting month progressed 
systematically with latitude. At some stations large differences were '* 
found in the same row in adjacent columns, so the procedure was repeated 
using averages of adjacent months. Only at Churchill and V^ite Sands did 
the latter e ’rcise result in larger maximum differences. 

Table 3 gives the differences of the mean magnitude of the 
surface temperature anomalies betvjeen the year groups following positive 
and negative zonal wind anomalies at 35-45 km. For brevity, only the column 
which contains the largest difference found is presented. The average 
values for the tv;o groups of months are presented in Figure 1 combined 
for all five stations. Testing against the null hypothesis indicates 
that the difference between the t\;o curves in Figure 1 is significant at 
the 57, level for months five and six. 

Grouping the months according to the sign of the anomaly in 
stratospheric vjind speed is admittedly tenuous, particularily when the 
difference betv;een a positive and negative anomaly may ho less than 1 m.s 
(e.g.. Table 4). However, this procedure yields results which are consis- 
tent auiong all stations used so is reported here in detail. In the absence 
of a jihysical e::plariati on tor the results, it woul * be more intuitively 
satisfying if they pertained to extreme anomalies, so the entire procedure 
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vas ropeaLed (after the fact) using only the months in which the extreme 
anomalies of zonal vind speed occurred. However, that exercise did not 
yield results having persistence from station to station and is not 
reproduced here. 

In an effort to stratify the surface temperature anomalies 
according to algebraic sign, tlic monthly means and anomalies of strato- 
spheric T and V were examined, along with the standard deviations and 
time and height derivatives of T, U, and V. None of those stratospheric 
parameters could successfully be employed to stratify the surface tempera- 
ture knomalics according to sign. 

Large surface temperature anomalies are usually associated with 
a persistent circulation feature, such as a blocking ridge (e.g.. Green, 

1969) with the algebraic sign of the anomaly determined by whether a 
station is upstream or dow’nstream. of the ridge. This suggests that the 
relationship reported above may be useful for predicting the occurrence 
of a" persistent circulation feature, but that single-station stratospheric 
parameters alone are not sufficient to determine at which longitude the 
ridge axis (or a similar phase-determining characteristic) will be found. 
Thus, a com.plete analysis of this relation should concentrate on the develop- 
ment of blocking higlis and other quasi-permanent circulation features, 
and should trace their origin to a particular pattern in the stratospheric 
circulation. Interlevel wave analysis technirpies are well-suited for this 
purpose and could be applied to the horizontal fields of rav;insonde, MRI'I, 
and satellite data taken from available maps. The negative zonal w’ind 
speed anomaly at a given station during the key month is related to a 
shift in phase and am.plitudc of the predominant standing wave features in 
the stratosphere. The amplitude and phase anomalies of zonal harmonics 
in the stratosphere during the key month at each station can be traced 
forw’ard in time (at the level whicli has highest coherence v;ith the 
preceding month) to the fifth and sixth months after the key month. The 
resulting mean time-height progression will cMii])i rically demonstrate the 
link between tlic zonal wind speed anomaly and the circulation features at 
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a later tiir.c at a particular level. Physical understanding of the chain 
of events leading to the tropospheric temperature anomalies will probably 
bo gained most readily by adapting one of the many numerical models avail- 
able in the literature and treating this as an initial value problem. By 
following the model output forward for five and six months the basic 
cause for blocking, or other quasi-permanent circulation features, to 
develop preferentially in one case may be detected. 

C. Correlation Studies of and Surface Anomalies 

Correlation coefficients betv/cen stratospheric anomalies of 
T, U, and V and surface temperature and precipitation anomalies were 
computed using monthly data for each of the individual twelve months, 
with the months grouped by season, and for the entire tim.e period. MilN 
data were used at nine individual height intervals, 20-60 km, and averaged 
in three layers (20-32, 34-46, 48-60 km). Correlation coefficients v;ere 
computed from lags zero through eleven months. Also, the analysis w?.s 
repeated using seasonal data (with four three-montli seasons \;here sums'cr 
is defined as June, July, and August) and correlation coefficients were 
computed from lags zero through seven seasons. Tlic following chart sum- 
marizes these correlation studies: 


S tations 


7 

Height Levels 

X 

12 

>1RN Parameters 

X 

3 

Surface Parameters 

X 

2 

Months /Seasons 


' 12 

X 

4 

Lags 


“ 12 

X 

8 


Total Correlations = 161,280 

Each value was statistically tested and tliosc which passed tlic 37o signifi- 
cance level were flagged. Tlie values were printed in a variety of fonr.ats 


to cxpidite comparison of levels or stations. Large correlation coeffi- 
cients were occasionally found, but no more often than expected by chance, 
and in most cases they v;ere not reproduced at another station nor did 
they show a pattern of progression with time or height. 

At Kennedy and Antigua the correlation coefficients between 
temperature at the surface and stratospheric T and u are persistently 
large (.Table 5) during the warm half-year (April-September) . However, 
it is well-known that surface temperature anomalies are persistent at 
tropical latitudes and, because of the quasi-biennial oscillation, 
stratospheric temperature and zonal wind anomalies also are persistent 
at low latitudes. Thus, while this result is consistent v;ith knowTi 
processes its usefulness does not seem significant. 

van Loon and Jenne (1975) report that persistence between seasons' 
is not generally a useful tool for forecasting surface temperature anom- 
alies. Ho\:evcr, as monthly data often show greater persistence than 
seasonal data, the above correlation study vjas repeated v?ith monthly data 
using the residuals of surface temperature and precipitation anomalies 
after autoregression at each lag. This procedure did not improve the 
results . 

IV. SUIlM/dlY 


In the present pilot study relationships between stratospheric 
anomalies or events and later surface temperature and precipitation 
anomalies at the same station have been sought. This single-station 
analysis has identified two relationships which may be useful for fore- 
casting purposes, and which may have intrinsic value for understanding 
the coupling betv;een the stratosphere and troposphere. The anomaly of 
precipitation received in October at stations north of 35 N is found to 
be related to the date the zonal wind changes direction at that station. 
Also, the magnitude of surface temperature anomalies five and six months 
following a key month are found to depend upon the stratospheric zonal 


vind speed during the key month at each station. IJotJi relations are 
statistically significant at the 5% level, but before either can be applied 
with confidence they must be analyzed further. The analysis of horizontal 
fields v:as beyond the scope of this study, but must be includeu in future 
studies designed to diagnose the preceding relationships. It seems likely 
that the anomalies used here reflect patterns in the large-scale circula- 
tion, both in the stratosphere and in the troposphere, and that under- 
standing of tropospheric-stratospheric coupling will eventually be in 
terms of the large-scale circulation. 

A further reason for analyzing horizontal and vertical fields, 
rather than isolated variables at a single station, is to take account 
of dynamic coupling processes. For example, several v.Titers have attemp- 
ted to link varying seasonal climatic regimes at mid-latitudes to the 
tropical quasi-biennial oscillation (most recently, Folland, 1975), and 
others (e.g., Winstanley, 1973) have tried to link low-latitude precipi- 
tation to the strength of the circumpolar vortex. Vertical coupling 
between atmospheric layers should be analyzed because the large-scale 
changes in upper level (20-80 km) circulation r..ay also influence the 
lower atmosphere just as it is known that the reverse is true, bugayeva, 
et al. (1975), report a slight increase in the success at the 500 mb level 
of a numerical prediction scheme when stratospheric forecasts arc included. 
The radius of correlaticn of meteorological variables generally increases 
with height, and a stratospheric link may be very important for under- 
standing atm.ospheric tcleconncctions and other seasonal climatic ciiangcs . 
Thus, inter level hemispheric or global patterns must be examined as a 
whole . 
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T\BLE 2. Anotnalies of the fall reversal at 30- km, October precipitation anomalies, and 200 mb 
heights. The reversal is defined as the last date the zonal wind’ became westerly 
after having been easterly for at least 5 consecutive days. (Negative anomaly means 
reversal v;as earlier than normal.) 
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TABLE 3. Differences of the mean magnitude of the anomaly in surface 
teinperature between years following a negative anomaly and 
years following a positive anomaly in zonal wind speed at 
35-A5 km. The selection of starting raonths is described in 
the text. Units are tenths of °C. The largest value in each 
column is underlined. 
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A. 

Mean monthly zonal wind speed anomali«'s at 35 

-A 5 km during 




the key 

month at Mills 

stations . 
in ms'^. 

Key months are defined in 




the text 

. Speeds are 

M is missing. 
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STATION 

MONTH 


196A 

1965 1966 

1967 1968 

1969 1970 

1971 1972 1973 

GREELY 

SEP 


-1.2 

5.2 1.9 

1.6 0.9 

-2. A -1.3 

M 1.7 -6.5 

CHURCHILL 

JUL/AUG 


1.0 

-0.5 0.6 

-0.8 0.1 

-O.A 2.0 

-0.3 -l.A -0.2 

WALLOPS 

JUN 


2.5 

0.6 2.5 

3.0 -1.1 

-2.5 1.6 

-5.0 -1.6 0.1 

FT’. MUGU 




1.8 

0.7 0.9 

2.1 -0.7 

-l.A 0.9 

-A. 9 -1.9 2.1 

WHITE SANDh 

IL 

\Y/JUN 

0.1 

-0.3 -0.3 

5.8 -1.3 

A.O -0.7 

-A.O -3.5 1.3 






TABLE 5. Correlation coefficients of anomalies in surface temperature with stratospheric (a) U 
and (b) T, based on monthly data for April - September. Lags, in months, are for 
stratospheric data with respect to surface data. The number of data pairs (N) is given 
only if a value exceeds the 57o significance level. Correlation coefficients are times 
one hundred. 
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Figure 1. Mean raagnitude of anoTnalics of surface temperature 

at five MRN stations up to eight r.ionths follov;ing the 
starting month (see tent). Dashed curve is for 
positive anomaly in zonal wind speed at 35-45 km 
the starting m.onth, solid curve is for negative 
anomaly in zonal v;ind speed. 



